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Abstract 

In-medium nucleon electromagnetic form factors are calculated in the quark 
meson coupling model. The form factors are typically found to be suppressed 
as the density increases. For example, at normal nuclear density and ~ 
0.3 GeV^, the nucleon electric form factors are reduced by approximately 8% 
while the magnetic form factors are reduced by only 1 — 2%. These variations 
are consistent with current experimental limits but should be tested by more 

precise experiments in the near future. 
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There is now considerable evidence that nucleons bound in an atomic nucleus experience 
very strong, effective scalar and vector helds M- It is a fundamental issue in nuclear 
physics to understand whether these strong helds alter the internal structure of the nucleon 
to a signihcant extent. There are already strong, experimental constraints on the extent 
of such changes [§|-§|. Furthermore, new quasi-elastic (e, e') and {e,e'N) experiments at 
facilities such as Mainz, MIT-Bates and especially TJNAF should improve these constraints 
in the near future. Our aim here is to examine the changes in the electromagnetic form 
factors of bound nucleons within one particular quark model of nuclear structure, the quark 
meson coupling (QMC) model In particular, we show that the model is consistent 

with existing constraints but that the predicted changes should be testable with only a 
modest increase in precision. 

As in Quantum Hadrodynamics (QHD) [|l|, the QMC model describes the properties of 
nuclear systems using effective scalar (cr) and vector (a;) helds. However, unlike QHD where 
the nucleon is a point particle, the nucleon in this model has substructure. The a and uj 
helds are coupled directly to the quarks within the nucleons, rather than to the nucleons 
themselves. As a result, the internal structure of the bound “nucleon” is modihed by the 
medium with respect to the free case. The relative simplicity of the QMC model at the 
hadronic level, together with its phenomenological successes, makes it suitable for many 
applications in nuclear physics ||I0 HI^] . 

The structure of the nucleon in the QMC model is described by the MIT bag [|T^. (How¬ 
ever, there are also versions which use the color dielectric model or a relativistic oscillator 
|T^ .) The small mass of the quark in these models implies that the lower component of the 
quark wave function responds rapidly to the change of its environment {a held), with a 
consequent decrease in the scalar baryon density. As the scalar baryon density itself is the 
source of the a held, this provides a mechanism for the saturation of nuclear matter where 
the quark structure plays a vital role. 

In this article we use the self-consistent solutions of the QMC model to calculate the 
density-dependent electromagnetic form factors of the nucleon. For simplicity, possible oh- 


2 




shell effects of the nucleon in medium |jT^ are ignored in the present treatment. The nucleon 
bag contains three independent quarks. With the static spherical cavity approximation to 
the MIT bag |^, the center-of-mass motion has to be corrected. We use the Peierls- 
Thouless projection method together with a Lorentz contraction for the nucleon internal 
wave function. This technique for implementing the center-of-mass and recoil corrections 
has been quite successful in the case of the electromagnetic form factors for free nucleons 


For an on-shell nucleon, the electric {Ge) and magnetic {Gm) form factors can be con¬ 
veniently defined in the Breit frame by 


(NAl)\Aom(-p ^xIx.Ge(Q\ ( 1 ) 

(A'u|)ldo)|A',(-|)> (2) 

where Xs and yj/ are Pauli spinors for the initial and final nucleons respectively, ^is the three 
momentum transfer, and = —Q^ = —The major advantage of the Breit frame is that 
Ge and Gm are explicitly decoupled and can be determined by the time and space compo¬ 
nents of the electromagnetic current (J^), respectively. Note that, in the above definitions 
[Eqs. (0) and ], both the initial and final states are physical states which incorporate 
meson clouds. 

In the QMC model, the quarks typically move much faster than the nucleon in the 
medium (the Fermi motion). Thus it is reasonable to assume that they always have suffi¬ 
cient time to adjust their motion so that they stay in the lowest energy state |^. In the 
instantaneous rest frame of a nucleon in infinite nuclear matter, the Lagrangian for the 
quark-meson coupling model is 


Cg = - mg)q9v 


BOv + glqqa 




1 1 
2 2 1 -*- 2 2 

) 


( 3 ) 


where niq is the current quark mass, B is the bag constant, gl and g^ denote the quark- 
meson coupling constants, and 6y is a step function which is one inside the bag volume and 
vanishes outside. In mean field approximation for the meson field, the normalized solution 
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for the lowest state of the quark is given by |T3 

/ 


where 


v47r 




jo{uJor/R) 
yi(3q(T-fii^ujQr/R) j 


9{R-r) Xm, 


6q -|- QijjUjR, 


Pg = 


fig — m*R 

\^q + m*qR' 


iVg = 2R^j^{uJo)[Vlq{flq - 1) + m*R/2]/ujl, 


(4) 

(5) 

( 6 ) 


with ilq = ^Uq + {m*Ry, m* = rUg — g^a, R the bag radius, and Xm the quark Pauli spinor. 
The frequency (cuq) of this lowest mode is determined by the boundary condition at the bag 
surface, 


jo(a;o) = pgip^o)- 


(7) 


The form of the quark wave function in Eq. is almost identical to the free case. However 
the parameters in the solution, Eq. (^, have been modihed by the medium. The mean 
values of the scalar held a and vector held uJ are self-consistently determined by the coupled 
equations of motion for the system In particular, a is given by the thermodynamic 

condition 


dE. 


total 


da 


= 0 , 


( 8 ) 


R,p 


where p is the baryon density and the total energy per nucleon for symmetric nuclear matter 
is given by 
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T^total 


(27r)^p 




and uj is determined through baryon number conservation to be 


_ ^9l 

uj = —p. 
m,, 


(9) 


( 10 ) 


The ehective mass of the nucleon in the medium can be expressed as 
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( 11 ) 


* _ ^0 I 

-^ + 2^-^ -O) 

where the second term on the right hand side of the above eqnation parameterizes the snm 
of the center of mass motion and glnon corrections. In the present work, the valnes of Zq 
and B are assnmed to be density-independent (the former for reasons explained in Ref. |P). 
They are determined by reqniring that the free nncleon mass be given by mAr = 939 MeV 
and by the eqnilibrinm condition in free space for a given bag radius, 


dm 


N 


dR 


= 0 . 


( 12 ) 


p=0 


For hnite nuclear density, the in-medium bag radius R can be obtained by solving Eq. (|T^ 
at hnite p - it typically decreases by a few percent at nuclear matter density (po)- 
The electromagnetic current of the quark is simply 


3q{x) = E<5/eg/(a:)7^g/(x), (13) 

/ 

where qf{x) is the quark held operator for the havor / and Q/ is its charge in units of e. The 
momentum eigenstate of a baryon is constructed by the Peierls-Thouless projection method 


0,0 


TpT(Xi,f2,^3;P) = iVpTe*^'^"'" g(Xi - fc.m.)g(x2 -Xc,m.)g(x3 -Xc.m.), (14) 


where A^pp is a normalization constant, p the total momentum of the baryon, and Xc.m. = 
(xi -f a ;2 -1- T3)/3 is the center of mass of the baryon (we assume equal mass quarks here). 

Using Eqs. ([I3| ) and (|^, the nucleon electromagnetic form factors for the proton’s quark 
core can be easily calculated by 

Ge{Q^) = J d^rjo{Qr)pg{r)K{r)/DpT, (15) 

Gm(Q^) = (2mM/Q) Jd^rji{Qr)(3gjo{LJor/R)ji{uJor/R)K{r)/DpT, (16) 

Dpt = jd\pg{r)K{r), (17) 

where Dpp is the normalization factor, Pg{r) = jQ^uor/R) + Pgjf{uJor/R), and K{r) = 

Jd?x pg{x)pg{—x — r) is the recoil function to account for the correlation of the two spectator 
quarks. 
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Apart from the center-of-mass correction, it is also vital to include Lorentz contraction 


of the bag for the form factors at moderate momentum transfer |]T6| , p!8| . In the prefered 
Breit frame, the photon-quark interaction can be reasonably treated as instantaneous. The 
hnal form of the form factors can be obtained through a simple rescaling, i.e.. 


Ge{Q'^) — 


(18) 

Gm{Q‘^) = 


(19) 


where E* = + Q'^/4: and are the form factors calculated with the static 

spherical bag wave function [Eqs. 


and ([T6|) ]. The scaling factor in the argument 
arises from the coordinate transformation of the struck quark and the factor in the front, 
{m*pf/E*Y, comes from the reduction of the integral measure of two spectator quarks in the 
Breit frame [ITB . 


As is well-known, a realistic picture of the nucleon should include the surrounding meson 
cloud. Following the cloudy bag model (CBM) we limit our consideration on the 

meson cloud correction to the most important component, namely the pion cloud. As in 
free space, the pion held is a Goldstone boson and acts to restore the chiral symmetry. The 
Lagrangian related to the pion held and its interaction, within the pseudoscalar quark-pion 
coupling scheme, is 

.n 1 




( 20 ) 


where 5s is a surface delta function of the bag, the pion mass and the pion decay 
constant. The electromagnetic current of the pion is 


jl^{x) = —ie[7r^{x)d^7r{x) — TT{x)d^Ti\x)], 


( 21 ) 


where n^x) = -^[xi{x) -|- i'K 2 {x)] either destroys a negatively charged pion or creates a 
positively charged one. As long as the bag radius is above 0.7 fm, the pion held is relatively 
weak and can be treated perturbatively. A physical baryon state can be expressed as m 


|2l) = 7^[1 + (m^ -H„- |2l„), 


( 22 ) 
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where A is a projection operator which projects out all the components of |A) with at least 
one pion, Hj is the interaction Hamiltonian which describes the process of emission and 
absorption of pions. The matrix elements of Hj between the bare baryon states and their 
properties are [p0| 


v^fik) ^ {Ao\Hj\7r^ik)Bo) 
w^fik) ^ {Ao7r^ik)\Hj\Bo) 


ifo 


AB 


u{kR) 


m, [2wt(2ir)3]V2 


T.CZ 


msA(^* (i* 

SsISa 


..BA 

Aj 


(k) 


= —V, 


AB 

Oj 


(k) = v^f(-k), 


(23) 

(24) 


where the pion has momentum k and isospin projection j, fA^ is the reduced matrix element 


for the ttBq —>• Aq transition vertex, u{kR) = 3ji{kR)/kR, Uk = y k'^ + and Sm and tn 
are spherical unit vectors for spin and isospin, respectively. The ttNN form factor [u{kR)] is 
fully determined by the model itself and depends only on the bag radius. The bare baryon 
probability in the physical baryon state, is given by 


^2 = 


1 + E 


7, 


2 

0 \ 


B 




1 

P 


dk k^vA{kR) 


n -1 


127 r 2 Jo u!k{mA — tub — oJkA 


(25) 


where P denotes a principal value integral. 

Up to one pion loop, the corrections to the electromagnetic form factors arising from 
the pion cloud are described by the following two processes shown in Fig. 1(b) and 1(c). 
The detailed expressions for their contributions can be found in Ref. [|T^] with the following 
substitutions —>■ m*, m^, and Jt^ab JIab- 

In principle, the existence of the tt and A inside the nuclear medium will also lead to some 
modihcation of their properties. Since the pion is well approximated as a Goldstone boson, 
the explicit chiral symmetry breaking is small in free space, and it should be somewhat 
smaller in nuclear medium ||^. While the pion mass would be slightly smaller in the 
medium, because the pion held has little effect on the form factors (other than Geii), we use 
m* = tUtt. As the A is treated on the same footing as the nucleon in the CBM, its mass 
should vary in a similar manner as the nucleon. Thus we assume that the in-medium and 
free space N — A mass splitting are approximately equal, i.e., m\ — ~ — uin- The 

physical ttAB coupling constant is obtained by ~ There are uncertain 
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corrections on the bare conpling constant , snch as the nonzero qnark mass and the 
correction for spnrions center of mass motion. Therefore, we nse the renormalized conpling 
constant in onr calculation, ~ 3.03, which corresponds to the usual tiNN coupling 
constant, — 0.081. In the medium, the nNN coupling constant might be expected to 

decrease slightly due to the enhancement of the lower component of the quark wave function, 
but we shall ignore this density dependence in a hrst treatment and use f*^j^ — fnNN- 

Fig. 2 summarises the nucleon electromagnetic form factors in free space with the 
bag radius R = 1 fm. The ‘dipole’ refers to the standard dipole £t, Fd{Q‘^) = 1/(1 + 
Q^/0.71 GeV^)^. The corrections for the center-of-mass motion and Lorentz contraction 
lead to signihcantly better agreement with data than was obtained in the original, static 
CBM calculations. For a detailed comparison with data we refer to Ref. |^. The purpose 
of Fig. 2 here is simply to show that our model produces realistic form factors in free space. 

Our main results, namely the density dependence of the form factors in matter, relative 
to those in free space, are shown in Fig. 3. The charge form factors are much more sensitive 
to the nuclear medium density than the magnetic ones. The latter are nearly one order of 
magnitude less sensitive. Increasing density obviously suppresses the electromagnetic form 
factors for small Q^. For a hxed (less than 0.3 GeV^), the form factors decrease almost 
linearly with respect to the nuclear density, p . At ~ 0.3 GeV^, the proton and neutron 
charge form factors are reduced by roughly 5% and 6% for p = O.Spo, and 8% for the normal 
nuclear density, po; similarly, the proton and neutron magnetic form factors are 1% and 
0.6% smaller for p = O.Spo, and 1.5% and 0.9% for the normal nuclear density. 

The best experimental constraints on the changes in these form factors come from the 
analysis of p-scaling data. For example, in Fe the nucleon root-mean-square radius cannot 
vary by more than 3% 0. However, in the kinematic range covered by this analysis, the 
eN cross section is predominantly magnetic, so this limit applies essentially to Gm- (As the 
electric and magnetic form factors contribute typically in the ratio 1:3 the corresponding 
limit on Ge would be nearer 10%.) For the QMG model considered here, the calculated 
increase in the root-mean-square radius of the magnetic form factors is less than 0.8% at po. 
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For the electric form factors the best experimental limit seems to come from the Coulomb 
sum-rule, where a variation bigger than 4% would be excluded |p]. This is similar in size to 
the variations calculated here (e.g., 5.5% for Gep at po) and not sufficient to reject them. 

In conclusion, we have calculated the density-dependent electromagnetic form factors 
of a bound “nucleon” within the QMC model. The differential cross section for the eN 
scattering can be easily constructed once the kinematics is selected. This work represents, 
to our knowledge, the hrst quantitative, density-dependent study within a quark model which 
reproduces the saturation energy, density and compressibility of nuclear matter. While the 
deviations of the form factors from their free values are within current experimental limits, 
they provide a strong motivation for forthcoming quasielastic electron-nucleus scattering 
experiments. 

We would like to thank I. Sick for helpful correspondence concerning the constraints from 
p-scaling. This work was supported by the Australian Research Council. A.W.T and K.S. 
acknowledge support from the Japan Society for the Promotion of Science. 
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FIGURES 


FIG. 1. 
pion loop). 



B 

(a) (b) (c) 


Diagrams illustrating the various contributions included in this calculation (up to 
The intermediate baryons B and C are restricted to the N and A. 


one 
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FIG. 3. The nucleon electromagnetic form factors in the nuclear medium (relative to those 
in free space case). The free space bag radius is 1 fm and the density is quoted in units of the 
saturation density of symmetric nuclear matter (po = 0.15 fm“^). 
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